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Abstract 

Two semi-analytical approaches to solve the problem of light scattering on nanowire antennas are developed 
and compared. The derivation is based on the exact solution of the plane wave scattering problem in case of 
an infinite cylinder. The original three-dimensional problem is reduced in two alternative ways to a simple 
one-dimensional integral equation, which can be solved numerically by a method of moments approach. 
Scattering cross sections of gold nanowire antennas with different lengths and aspect ratios are analysed for 
the optical and near-infrared spectral range. Comparison of the proposed semi-analytical methods with the 
numerically rigorous discrete dipole approximation method demonstrates good agreement as well as superior 
numerical performance. 
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1. Introduction 

Since decades antennas are used in everyday de- 
vices in the radio and microwave spectral range as a 
bridge between propagating radiation and localized 
fields. In this spectral range semi-analytical models 
exist [l|, offering insight into interaction processes 
and guiding engineers in antenna design. In the 
same time, for practical antenna design and opti- 
mization well established numerical tools are typi- 
cally used • Shifting antenna resonances towards 
optical spectral range brings new challenges both 
from the fabrication and the theoretical perspec- 
tives. At optical frequencies metal can no longer 
be treated as perfect electric conductor and dimen- 
sions of the antenna might be as small as several 
tens of nanometers Q. Recent progresses in nan- 
otechnology have enabled the fabrication of optical 
antennas (nano-antennas) Q , Q and opened many 
exciting possibilities towards nano-antenna appli- 
cations. For example, it has been recently demon- 
strated that nano-antennas can enhance |5J and di- 
rect the emission of single molecules Q and that 
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they can play a key role in sensing application 0]. 
Great potential in improving the efficiency of solar- 
cells should also be mentioned [8|. 

Design and optimization of optical antennas 
are mainly done using general numerical Maxwell 
solvers which demand huge computational re- 
sources. Therefore accurate analytical and semi- 
analytical models predicting characteristics and 
performance of nano-antenna are of great impor- 
tance. There are just very few exact analytical so- 
lutions available. Light scattering on metal spheres 
[Io| . infinite long cylinders and spheroids 
can be derived in closed analytical form. For the 
light scattering problem involving a finite length 
nanowire, a Pocklington-like one-dimensional (ID) 
integral equation can be introduced using equiva- 
lent surface impedance method |12]. In this pa- 
per, we propose a further development of this ap- 
proach, which provides better accuracy especially in 
the case of nano-antennas with small aspect ratio. 

The paper is organized as follows. In section 2 
the problem of light scattering on a thin perfectly 
conducting wire is reviewed. Pocklington's integral 
equation is introduced and extended to the case of 
a nanowire of finite conductivity. We demonstrate 
how one can use the knowledge of the exact solution 



Preprint submitted to Elsevier 



January 15, 2013 



7\ 



E 



Figure 1: Definition of the geometrical parameters, radius a 
and length I, of the scattering cylinder as well as the chosen 
body centered coordinate system. Additionally the incident 
angle £ and the polarization basis vector Eljj are depicted in 
the incident plane. 



of the problem of the plane wave scattering on an 
infinite cylinder in order to imp rove accuracy of the 
surface impedance method [12J . In Sections 3 and 4 
a new numerical method to solve the resulting one- 
dimensional (ID) integral equation is introduced. 
The method involves a method of moments (MoM) 
like discretization scheme, but does not require any 
specific boundary conditions to be imposed at the 
nanowire ends. In section 5 numerical calculations 
of scattering cross-sections for plane wave scatter- 
ing on gold nanowires with varying geometries is 
presented and compared with numerically rigorous 
discrete dipole approximation (DDA) calculations 
[l^ . Section 6 summarizes the paper. 



2. Pocklington like equation 

In infinite free space the electric field E(r) gener- 
ated by a time harmonic current density distribu- 
tion j(r) (time dependence e~ lu3t ) enclosed in the 
finite volume V is given by (l4l | 



E(r) = E""-(r) H- wn (¥ + pV®V 



(r,r')j(r')dV (1) 
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where I denotes the three-dimensional unit ten- 
sor, ® is the tensor product defined by (a® b)^ — 
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is the scalar Green's function and k — — the free 

c 

space wave number. E mc is an electric field due 
to sources not contained in V. In scattering prob- 
lems the driving current density is not controlled 
from the outside, but instead induced by E mc which 
plays the role of an incident field in this case. Spa- 
tial derivatives of the integral in fTJ) exist as long as 
the source current density j(r) satisfies the Holder 
condition (Hi, |j (r) - j(r')| < fc|r - r'| Q with k > 
and < a < 1 for every pair r, r' 6 V. 

In what follows we consider light scattering on a 
cylinder (wire) with length I and radius a (Fig. [TJ) . 
The geometry of the problem including a body cen- 
tered coordinate system with z-axis parallel to the 
cylinder axis is depicted in figure [1] The wave vec- 
tor k enclosing the incident angle £ with the positive 
z-axis lies in the sz-plane (incident plane). To ex- 
cite longitudinal resonances only the projection of 
the incident plane wave electric field on the incident 
plane Eg = Eq (sin£z — cos£x) have to be taken 
into account. The perpendicular polarization com- 
ponent can be safely ignored. Under the assump- 
tion that the cylinder diameter is much smaller than 
the free space wavelength, i.e. ka <C 1, the incident 
electric field interacting with the cylinder can be 
viewed function depending only on z 
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First we briefly review the derivation of Pock- 
lington's equation for scattering on a thin perfectly 
conducting wire [lj. In this case in cylindrical co- 
ordinates {p, (j), z} the induced current density j(r) 
has only a z-component, shows no (^-dependence 
and exists solely on the antenna interface. Then 
the induced current I(z) is related to the current 
density j (r) as 



j(r) = zl(z 
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so that 



I(z) = 
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d(j) / pdpj(r). 
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Using ((U) in ([lj) and performing the //-integration 
one yields for the z-component of the electric 
field on the cylinder surface the following integro- 
differential equation 



E z (a,z) = E l z nc (z)+i 
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where g a (<t>', z — z') — g(r, r') with r = (a, 0, z) and 
r' = (a, <f)' , z') in cylindrical coordinates. While the 
antenna is perfectly conducting, the ^-component 
of the total field at the antenna surface has to van- 
ish, E z (a,z) = 0. Applying this boundary condi- 
tion to Eq. © one can derive Pocklington's integro- 
differential equation in the form 



i_ j? 2 

fc 2 dz 2 



2tt 



dz' d<t>'g a (ct>',z-z')I(z'), (7) 
Jo 



which provides a general solution of the scattering 
problem Here the fact, that the antenna is per- 
fectly conducting, is taken into account twice, first 
by assuming a special form of the induced current 
(HJ) and second by enforcing the boundary condition 
E z (a,z) = 0. 

A typical nano-antenna at optical frequency 
range possesses high, but finite conductivity. In 
this case the surface current approximation 
is still applicable, while the boundary condition 
E z (a,z) = is generally not. In order to use a 
Pocklington's like equation at this frequency range, 
one needs a relationship between the field E z {a 1 z) 
and the current I(z) at the antenna interface. 

Assuming a long nanowire, i.e. / 3> a, it is rea- 
sonable to expect that the internal electric field is 
separable similar to the solution of the equivalent 
problem involving an infinite cylinder [10]. More- 
over, it can be shown that as long as the wire is 
electrically thin, i.e. ka <C 1, and |e r | 3> 1 at least 
for the case of inclined incidents (£ ^ §) the p- 
and ^-components of the internal field are negligi- 
ble compared to the z-component of the field and 
the z-component shows no (^-dependence. Assum- 
ing these conditions arc fulfilled the internal electric 
field can be approximated by 
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where Jo (k p p) denotes the Bessel function of the 
first kind, k p — k\Je r — cos 2 £ with relative per- 
mittivity e r of the wire and f{z) an unknown z- 
dependent function giving the amplitude of the in- 
ternal field along the wire. A connection between 
the induced current density and the internal electric 
field is given by means of the volume equivalence 
theorem [lj by 



with Ae r = € r — 1. Combining ((5J) and © the total 
current through the wire can be calculated 
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Further comparing results of the integration in (JTI 
with ([5]) one can derive the following relation be- 
tween the electric field and the total current at the 
wire interface 



E z (a,z) = Z s I(z) 
with the surface impedance 

J {k p a)k p 



Z s = i 



2naujeoAe r Ji(kpa) 



(11) 



(12) 



Using relation (TTTJ) the following Pocklington like 
integro-differential equation for induced total cur- 
rent, the surface impedance (SI) integro-differential 
equation [12], can be obtained from equation © 



Z s I(z) = E l z nc (z) + i 



up 
2tt 



1 



1_^_ 

fc 2 dz 2 



dz' / d(j)' g a ((j)' , z 
Jo 



z')I(z'). (13) 



We propose a further improvement to the approx- 
imation (|13p by releasing the solely surface current 
ansatz ((4]). In order to do that, we assume that the 
induced current density on the right hand side of 
equation (TTJ) can be factorize similar to the internal 
field ©• In this way substituting ([S} in equation 
(TT]) both on the left hand side as boundary condi- 
tion as well as by using (0) on the right hand side 
to rewrite the induced current density one obtains 
a self-consistent integro-differential equation for the 
unknown amplitude f(z) 
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where g(a,z;r') = g(r, r') with r = (a,0,z) 
in cylindrical coordinates. This volume current 
(VC) integro-differential equation takes into ac- 
count both appropriate boundary conditions at the 
wire interface and an appropriate volume current 
distribution inside the wire. 
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3. Discrete form of VC integral equation 

In order to solve numerically the integro- 
differential equations (TT51) and (fT4")l one has to im- 
pose additional boundary conditions at the nano- 
antenna edges. A common choice is to impose the 
total current I(z) to be equal to zero for z = ±1/2 
[l^ . For a solid wire with finite conductivity this 
choice is generally not justified, while the total cur- 
rent can be discontinuous at the wire edges 

To overcome the requirement of additional 
boundary condition one has to convert the integro- 
diffcrential equations into purely integral ones. To 
do that one has to bring the differential operator 
(l + k~ 2 d 2 ) inside the integral. This procedure re- 
sults in singularities of the order |r — r'| 3 which 
are generally not integrable over a volume. To 
treat these singularities we follow the regulariza- 
tion scheme proposed by Lee et al. in Ref. (l6| . 
The regularized VC equation (|14p reads: 

(l + Ae r L 33 )E z (a,z)=ET c (z)+ 

lrA< : I I G 33 (a,z;r , )E z (p , ,z')d 3 r'+ 
'v-v* 
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with E z (p, z) = f(z)Jo (k p p). Here V* denotes a fi- 
nite and arbitrary shaped principal volume V* con- 
taining the singular point r = (a,0,z), go{r,r') — 
linife^o ff( r J r ') is the static scalar Green's function, 
G 33 is the zz-element of the dyadic Green's function 

a 



G 33 (r,r')=(l + ^) 5 (r,r'), (16) 
and L 33 the zz-element of the source dyadic [3 
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The surface integration in (I17[) has to be performed 
over the surface dV* enclosing the principal volume 
V*, n denotes the outer surface normal. 

The main advantage of the regularization scheme 
[3 is ( i) that all singularities disappear and (ii) 
that the principal volume V* can be finite and ar- 
bitrary shaped. We choose a cylinder with length 
A, radius y and center at the singular point r = 



(a, 0, z) as the principal volume (Fig. [2]). We assume 
A to be small, such that both the Bessel function 
Jo {k p p') and the amplitude function f(z') are ap- 
proximately constant over V*. Taking that into 
account, writing the volume integrals in the reg- 
ularized VC equation (fT5"|) in cylinder coordinates 
{p, <f>, z} and collecting all terms containing f(z) on 
the left and all terms containing f(z' ^ z) on the 
right hand side, one obtains a one-dimensional in- 
tegral equation in the form 



f(z)T^Er(z)+ J dz'f(z>)£(\z-z'\)+ 
[ 2 dz'f(z')£(\z-z'\), (18) 
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and 
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In (|2"0|) the integration is performed over the princi- 
pal volume V*, while in (|2ip over the corresponding 
wire slice of thickness A centered at z = z' but with 
excluded principal volume V*. To calculate these 
two integrals a new coordinate system with its cen- 
ter at the singular point has been chosen Fig. @. 
In the new coordinate system the source dyadic (|17|) 
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can be explicitly written as 
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The radius vector p' and the integration ranges in 
(|2QI21l23p parametrically depend on the nanowire 
radius a and the wire slice thickness A and are given 

by 

' ' a \ / pcosi; 
/ \ psint/ 



(24) 
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Taking into account that the nanowire diameter 
is small in comparison to the wavelength one can 
further simplify the integrals in (|2"0)) and (f2"Tj) by 
expanding the Green's function G33 in power of k 
up to the linear term 

J_ ( ZR\-R 2 1 

4tt I, k 2 R 5 + 2R 3 
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(26) 

where R = r' — r. Using (|26[l one can analytically 
integrate equation (|20p over z and p and equation 
(|2"Tj) over z . Residual integration in (|2T)|) , (f2"Tj) and 
(|23p has to be performed numerically. The regu- 
larization scheme (|15p ensures, that numerical in- 
tegration converges. In this way T in (|T5|) can be 
efficiently calculated once for given radius a of the 
nanowire antenna, A and wavelength. 

Equation (fT5|) can be solved numerically by a 
method of moments (MoM) approach. Specifically 
we choose a point matching MoM scheme with 
pulse- function basis [17j| . That means we divide 
the nanowire into a set of slices with thickness A 
and label each slice with an index i. For reasonable 
thin slices the amplitude function f(z) on each of 
them can be approximated by its value at the cen- 
ter f(zi) — U- This leads to the n — ^ dimensional 
matrix equation as discretized version of (|18|) 
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Figure 2: Cross-section of the nanowire (big circle) and the 
principal volume V* (small circle) together with the defini- 
tion of the coordinate systems used to calculate integrals in 
equation (1181) . 



which can be easily inverted numerically to yield 
the discrete set of values f{zi). While the Green's 
function depends on the distance between two slices 
\z — z'\, only the first row M\j of the matrix M has 
to be calculated. All other elements can be filled 
using the rule My = Mj-ij-x- 

Using the Green's function expansion ([26)) in the 
calculation of the matrix elements My for j < =^ 
the integration over z' can be done analytically. 
Only the integration over the nanowire profile has 
to be done numerically. In the calculation of the 
matrix elements M\j for j > ^ the full dyadic 
Green's function have to be used and so the com- 
plete volume integral has to be done numerically. 
These integrations are done over r'-regions far from 
the singularity and they demonstrate good con- 
vergence. An additional performance improvement 
can be achieved by expanding the Bessel function 
Jo for small arguments as [18] 
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Having the amplitude function f(z) calculated from 
([2"Tf the z-component of the induced electric field 
is given by the ansatz ([5]). 

4. Discrete form of SI integral equation 

To solve the surface impedance integro- 
differential equation (|13p the regularization 
and discretization scheme presented in section [3] 
can be also applied. A principal volume is chosen 
as a full slice of the wire at position z with 
thickness A. Regularizing the SI equation ([T3")) . 
writing it in cylinder coordinates and collecting all 
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terms containing I(z) and I(z' 7^ z), one obtains a 
one-dimensional integral equation in the form 



i(z)r' 
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dz'I(z')£ s (\z- z'\)+ 
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and G3 3 (0', z — z') = G33(r, r') as well as g§((/)', z — 
z') = limfc^o §( r , r ') with r = {a, 0, 2;} and r' = 
{a, </>',z'} in cylinder coordinates. The source 
dyadic if 3 can be expressed in terms of the com- 
plete elliptic integral of the second kind E (to) [IH 
as 

1 8 / 16a 2 * 



1,33 7T 16a 2 + A 2 
where E (to) is defined by 
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In r s ((22) the z'-integrations can be performed 
analytically with the help of the near-field expan- 
sion of the Green's function (|2"6| . The residual 
(^'-integration has to be done numerically once for 
given nanowire radius a, A and wavelength. 

To solve (|3"0"|) numerically we use a point match- 
ing MoM scheme with pulse function basis [ItJ as 
in section [31 Again the nanowire should be divided 
into a set of slices of thickness A yielding the matrix 
equation 



(36) 
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Figure 3: Scattering cross-section of a gold nanowire (I = 
200ram, a = lOnm) under slanting incidence (§ = -j) cal- 
culated with different rigorous numerical methods as well 
as the proposed volume current and surface impedance one- 
dimensional integral equations. In the top panel (a) the first 
and in the bottom panel (b) the third resonance peak are 
shown. 



where 
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for i =/= j 
for i = j 
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and Ii — I(zi). Similar to the case of volume cur- 
rent integral equation considered in section [3] only 
the first row of the matrix rvt s have to be calcu- 
lated due to the symmetry of the Green's function. 



Further for slices with \z\ 



< 2a the near field 



approximation of the Green's function (|26[) can be 
used to perform the z'-integrations in (|37p analyt- 
ically. The remaining integrations as well as the 
matrix inversion have to be performed numerically. 
Good convergence of the integrals is ensured, while 
the integrand is evaluated far from the singular 
point. Numerical inversion of Q36p results in the 
current Ii — I(zi) at the discrete number of points 
Zj along the nanowire. Finally the z-component of 
the induced electric field is given via ([TTj) . 

5. Numerical results and discussion 

In this section the semi-analytical methods de- 
veloped in sections 3 and 4 are evaluated and com- 
pared with numerically rigorous methods. Plane 
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wave scattering on a gold nanowire is considered in 
the optical and near-infrared spectral range. The 
relative permittivity of gold in this spectral range 
can be modeled by a free electron Drude model 



e r (u) 



with parameters e c 
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and 7 = 7.59297 • 10 13 s 1 [9j. For comparison pur- 
poses we calculate the total scattering cross-section 



O scat 



,d 2 r |E S 

| E mc| 2 



(39) 



where the scattered electric field E s is generated 
by the induced current density j(r) given by equa- 
tion © in the case of the VC and by equation ^ 
in the case of the SI integral equation model, re- 
spectively. In the far-field zone the scattering field 
is purely transverse and one can use the far-field 
dyadic Green's function 



lim t?( 



iki r — r-r 



r, r 



4irr 



| I -r®f| (40) 

and integral relation 

E s (r)= lW/ i / d 3 r'H(r,r')j(r') (41) 

to calculate the scattering cross-section. For the VC 
integral equation model the scattering cross section 
can be written in terms of the amplitude function 

/(*) 
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while for SI integral equation model in terms of the 
total current I(z) 
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We have performed a systematic convergence test 
of the developed one-dimensional integral equa- 
tion based methods for gold nanowires of different 
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Figure 4: Scattering cross-sections of gold nanowire with 
fixed radius a = lOnm and different lengths. The full and 
dashed lines show the VC and SI integral equation results, 
respectively. The dotted line is the numerically rigorous ref- 
erence calculated using DDA method. 
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Frequency [THz] 



600 



Figure 5: Contours plot of the resonance frequencies of the 
nanowires with different lengths under slanting incidence 
(£ = J). Top panel radius is a = lOnm, bottom panel 
a = 20nm. The full and dashed lines show the VC and SI 
integral equation results, respectively. The dots are repre- 
sents the numerically rigorous DDA calculations. 
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lengths and radii. Spatial resolution of A = lnm 
(thickness of the individual slices in the wire dis- 
cretization) typically results in less than ITHz dis- 
crepancy with respect to the converged value of the 
scattering cross-section maximum calculated at res- 
olution as small as A = 0.05nm. In what follows 
spatial resolution of A = lnm has been used for all 
calculations based of the VC and SI integral equa- 
tion models. 

In figure [3] the first (top panel) and the third 
(bottom panel) resonance peaks of the cross-section 
spectra are shown for a gold nanowire (a = 
lOnm, I = 200nm) under slanting incidence (£ = 
j) with |E* nc = 1. Results of both different rigor- 
ous numerical methods and semi-analytical ID ones 
are compared. For the rigorous numerical calcula- 
tions we used (i) HFSS, a commercial finite-element 
frequency-domain Maxwell solver from ANSYS 
fl9^ . (ii) an in-house implementation of the finite 
difference time domain method (FDTD) [20| and 
(iii) ADDA, an open-source software package for 
calculating scattering parameters using the discrete 
dipole approximation (DDA) algorithm [21]. The 
space discretization in the shown DDA and FDTD 
calculations was set to lnm. Discrepancies among 
the scattering cross-section spectra calculated using 
different rigorous three-dimensional (3D) Maxwell 
solvers are comparable with discrepancies of these 
spectra with the spectrum calculated using VC in- 
tegral equation method. The accuracy of the SI 
integral equation method is slightly worse but still 
very reasonable. 

Most important are the differences in used com- 
putational resources and execution time. The cal- 
culation of one frequency point in VC (SI) inte- 
gral equation method requires approximately 2 (1) 
seconds on one core of a workstation using Mathe- 
matica [22|. In contrast DDA calculations requires 
around 8 minutes per frequency point on the same 
workstation. HFSS needs around 6 minutes per 
frequency point if the mesh is optimized for one 
frequency only and is reused without optimization 
for 30 other frequencies. With FDTD one gets the 
complete spectra in one run in around 250 minutes 
on one core. In conclusion the newly proposed ID 
semi-analytical methods provide a speed-up in exe- 
cution time close to 200 times compared to general 
3D Maxwell solvers. Additionally up to 100 times 
less RAM is required for the semi-analytical calcu- 
lations. 

In figure 0] the position of the scattering cross- 
section maxima are shown for gold nanowires with 



different aspect ratios. Radius is a = 10nm and 
length varies from I — 50nm (shown in the in- 
set) to I = 300nm. Normal incidence (£ = §) 
with |E mc | = 1 is considered. VC integral equa- 
tion method (full black line), SI integral equation 
method (dashed line) and DDA (dotted line) are 
compared. As it is expected the accuracy of both 
semi-analytical methods deteriorates with decreas- 
ing aspect ratio. However even for considerable 
low aspect ratio | (I = 50nm) the first resonance 
peak predicted using approximate methods repre- 
sents the rigorous numerical result with a good rel- 
ative accuracy (2% of the central frequency). 

In general VC integral equation method demon- 
strate better agreement with the DDA calcula- 
tions in comparison with the SI method (Fig. [3] 
and [5]) . This can be best seen in figure [SJ where 
the resonance frequencies of increasing order (from 
left to right) against the wire length are depicted 
for nanowires with radius a — lOnm (top panel) 
and a = 20nm (bottom panel). The VC integral 
equation method (solid line), SI integral equation 
method (dashed line) are presented. Dots repre- 
sent DDA results. The brighter (darker) curves are 
the resonances of even (odd) orders. All calcula- 
tions are done under slanting incidence (£ = ?). 
The better accuracy of the VC integral equation 
method can be systematical traced in figure [S] espe- 
cially for the structures with smaller aspect ratio. 

6. Conclusion 

Two alternative methods to solve the scatter- 
ing problem on optical nanowire antenna, the vol- 
ume current integral equation (VC-IE) method and 
the surface impedance integral equation (SI-IE) 
method are introduced. In order to reduce the 
general 3D volume integral equation describing the 
scattering problem to a simple semi-analytical ID 
integro-differential equation, both methods utilize 
solutions of the problem of plane wave scatter- 
ing on infinite cylinder. A regularization and dis- 
cretization scheme is proposed in order to trans- 
form integro-differential equations into solely inte- 
gral equation. This transformation enables to solve 
the original problem without necessity to impose 
additional boundary conditions at the nanowire 
edges. Numerical evaluation of the proposed meth- 
ods and their comparison with different numerically 
rigorous methods is presented for scattering cross- 
section calculations. Gold nanowires are analyzed 
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at optical and near-infrared spectral range. The in- 
troduced one-dimensional semi-analytical methods 
demonstrate good agreement and superior numer- 
ical performance in comparison with rigorous nu- 
merical methods. 
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